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constant parameter defined by equation (41), A = - — — 
optimization function, equation (37) 
iteration function, equation (38) 
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dimensionless heat transfer rate defined by equation (43), 
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ANALYSIS OF CONDENSATION ON A HORIZONTAL CYLINDER WITH 
UNKNOWN WALL TEMPERATURE AND COMPARISON WITH THE 
NUSSELT MODEL OF FILM CONDENSATION 


Parviz A. Bahraini 
Ames Research Center 


SUMMARY 


Theoretical analysis and numerical computations are performed to set forth a new model of film 
condensation on a horizontal cylinder. The model is more general than the well-known Nusselt 
model of film condensation and is designed to encompass all the essential features of the Nusselt 
model. It is shown that a single parameter, constructed explicitly and without specification of the 
cylinder wall temperature, determines the degree of departure from the Nusselt model, which 
assumes a known and uniform wall temperature. It is also shown that the Nusselt model is reached 
for very small, as well as very large, values of this parameter. In both limiting cases the cylinder wall 
temperature assumes a uniform distribution and the Nusselt model is approached. The maximum 
deviations between the two models is rather small for cases which are representative of cylinder 
dimensions, materials and conditions encountered in practice. 

INTRODUCTION 


Condensation of saturated vapors on solid surfaces under conditions encountered terrestrially and 
during space flight is of great importance in application of heat and mass transfer analyses. The early 
heat and mass transfer literature for condensing fluids contains many insightful analysis of the 
physics involved, as evident in reference 1. The well known Nusselt model of film condensation on 
a horizontal cylinder is perhaps the most important one, owing to its fundamental geometry. It 
assumes that a uniform wall temperature distribution in the cylinder prevails. Predictions of the heat 
transfer rates are based upon this temperature which is presumed to be known, (ref. 2). However, the 
cylinder wall temperature for this problem is generally not known, and one would rather deal with 
the bulk temperature of the cooling fluid inside the cylinder. The dynamics of the heat and mass 
transfer processes determine the wall temperature distribution and, therefore, the wall temperature 
can not be specified a priori, as assumed in the Nusselt model. 

The aim of the present analysis is to set forth a new and more complete model of film condensation 
on a horizontal cylinder, where no assumptions are made to explicitly specify the wall temperature. 
The new model will be compared with the Nusselt model to determine the extent of differences in 
the heat and mass transfer predictions for conditions most widely encountered in horizontal cylinder 
condensers. 



THE ANALYSIS 


The velocities in condensing fluids are generally very small such that, the inertial forces may be 
ignored in comparison to other forces, Eckert (ref. 3). Furthermore, the governing equations are 
distinctly boundary-layer in nature, Schlichting (ref. 4). A balance between the pressure, viscous and 
body forces within the condensing film prevails and the pressure within the condensing film is that 
of the quiescent saturated vapor away from the cylinder. The thickness of the condensing film is 
much smaller than the cylinder radius. Therefore, it is convenient to employ a curvalinear coordinate 
system in which the abscissa x lies on the cylinder surface and the ordinate y is perpendicular to and 
measured from the surface. 


With these assumptions, the Navier-Stokes equation of motion for this problem can be reduced to: 
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Denoting the distance from the uppermost point on the cylinder directly downwards by z, 
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one may write, 
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But, as noted earlier, the pressure within the condensing film is that of the quiescent saturated vapor 
away from the cylinder, 

dp 
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dx 


By combining equations (1) and (4), 


R— r = -(p-Pv)gsin<P 

dy 


and denoting the density difference by p , 
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P = P“Pv 
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equation (5) can be integrated: 
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( 7 ) 


02 . 0 

“ u = ^-sincpy + Cjy + c 2 

The no-slip boundary condition at the outer tube wall, u = 0 at y = 0, eliminates the second 
integration constant, c 2 . The shear force at the liquid-gas interface is taken to be negligible due to 
small viscosity of the vapor relative to that of the liquid, leading to evaluation of the first integration 
constant, Cj . 
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The mass flow rate in the condensed layer at any x is: 
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so that, 

dih = -2Md|5 3 sincpj (12) 

3fi 

is the incremental increase in the mass flow rate at any x. The heat transfer rate is simply, 

dQ = hfgdm (13) 


Using the series resistance concept for the fluid layer, the cylinder wall and internal heat transfer, 
and assuming that the internal heat transfer coefficient hj is an appropriate mean value, one may 
write, 
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By carrying out the differentiation in equation (12) and using equations (13) and (14) we get, 
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The boundary condition is due to the symmetrical film thickness at the top of the of the cylinder, i.e. 

(16) 


d5 _ 
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Now, by introducing a dimensionless film thickness A 
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and a dimensionless parameter O, 
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equation (15) gives, 
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Rearranging equation (19) results in a first order nonlinear ordinary differential equation, in a 
suitable form for Runge-Kutta numerical integration: 
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The boundary condition, equation (16) is transformed into, 
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dcp 


= 0,cp = 0 


( 21 ) 
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which, when used in the differential equation (19) gives, 

Tj- (22) 

1 + 

A 0 

Here, Aq denotes the value of A at cp = 0 and can be evaluated once the parameter Q is specified. 
At this point it is apparent from reference 2 that equation (20) is identical to the one obtained from 
the Nusselt model, when Q. — » 0. 


Recasting the Nusselt Model in Terms of the Variables of the Present Analysis 

In order to make a precise comparison of the predictions of the present model to those of the Nusselt 
Model, one must first recast the Nusselt model in terms of the variables of the present model. If a 
relationship between T b in the above analysis and T w in the Nusselt model could be obtained, one 
could compare the numerical solutions of equation (20) to the results of the Nusselt model. 


Applying the series-resistance concept again. 
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But, 


Q = m^lifg 

where, m n is the mass-flow rate as (p — » n, and by using equation (11), 
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A limit is imposed since, the mass flow crossing cp = k is zero and not a representative of the 
condensed mass flow rate. 

Equation (26) can be written in terms of A Nu and T w , by noting that in the present model as 
k w , hj^oo, T w — > T b , a constant and in effect, the Nusselt’ s model is approached. Therefore, T w 
replaces T b for A = A Nu . 
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Rearranging, 
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Using the definition for Q and introducing a dimensionless temperature % , 
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one obtains the following equation which will be utilized later in the analysis: 
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Writing equation (28) in terms of A Nu , % and (T sat - T^), one obtains: 
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Similarly for equation (26), one obtains: 
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where, in tandem with equation (30), 


A = lim(A 3 sin(p) (35) 

(j) — >71' ' 

We note that comparison of Q Nu and Q can now be made, once the value of Q is specified and the 
integration of equation (20) is carried out. 


Numerical Solution 

A Runge-Kutta numerical integration technique was used to integrate equation (20) from q> = 0 to 
(p = 179.99° . This method is designed to approximate the Taylor series method without requiring 
explicit definition or evaluations of derivatives beyond the first. The approximation is obtained at 
the expense of several evaluations of the function. Reference 5 describes the method in detail and 
provides a framework for a subroutine for solving initial value problems of ordinary differential 
equations. A modem modification to the classical Runge-Kutta technique has been employed to 
control the step size . Since only one solution value is required for calculation of the next, the 
method is self-starting. It requires six function evaluations per step. Four of these function values are 
combined with one set of coefficients to produce a fourth-order method, and all six values are 
combined with another set of coefficients to produce a fifth-order method. Comparison of the two 
values yields an error estimate which was used for step size control. The boundary condition at 
(p = 0 appears as an initial condition, where integration over a specified interval was performed in 

one or many subintervals. Specification of — = 0,® = 0 is not an essential part of the routine, 

d(p dA 

however, since the denominator of equation (20) contains a singularity at <p = 0, — was set to zero 

dcp 

at this point. Once Q. was specified, Aq was calculated and the integration was carried out. Local 
error of 10 5 and 10 8 yielded identical values of A up to the seventh significant digit. 


Calculation of A 0 

Equation (20) is a quatic equation and in general may possess multiple roots. The roots, however, 
may be real or complex conjugates. Real negative roots are not physical and complex roots are not 
within the scope of consideration of the present work. Here interest lies with the particular solutions 
(i.e., positive real roots) which enable comparison of the present model with the Nusselt model of 
film condensation. Inspection of equation (22) reveals that two real and two complex solutions can 
be obtained by considering its resolvent cubic equation, as described in reference 6. Furthermore, the 
real positive roots are equal. During attempts to integrate equation (20) for £2 > 0, it was discovered 
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that Aq had to be specified to a high degree of accuracy to avoid incompatibility between the initial 
value and the differential equation. A successive approximation scheme was devised for this purpose 
by letting. 


f(A 0 ) = At 
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where Ajj denotes the n-th approximation to A 0 , the function g(Ag) is the iteration equation and A 
is a constant. Then, 
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For g(Ao) = 0 one may choose Aq = 1 as a starting trial value, and convergence of the 

dAg 

successive approximations can be expected. 

4 A + 3 AQ + 1 = 0 


(39) 


Therefore, 


A = 


-1 
4 + 3Q 


(40) 


Specification of 12 

In order to specify values of 12 which are representative of cylinder dimensions, materials and 
conditions encountered in practice, a simple survey was conducted. Cylinder diameters between 
0.625 and 1.25 inches and wall thicknesses of 0.045 and 0.065 inches are most commonly used in 
terrestrial applications of horizontal tube condensers. The pipe materials are usually admiralty brass 
or stainless steel. It was necessary to obtain only approximate or orders of magnitude of various 
parameters. Standard correlations of turbulent pipe flow, Eckert (ref. 3), were used to evaluate the 
Nusselt number Nu, the cylinder internal heat-transfer coefficient hj and finally, the parameter, £2. 

For conditions of atmospheric pressure for saturated water vapor and Reynolds numbers, based on 
cylinder diameter, ranging between 10 4 and 10 5 , various calculations were carried out. A sample of 
the results is presented in Table 1 to highlight the expected range of values for Q. 
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Table 1. 


Re 

Nu 

hi W/K-m 2 

£2 

10,000 

104 

2.561 

3.05 

20,000 

185 

4,558 

1.73 

30,000 

260 

6,412 

1.24 

40,000 

332 

8,181 

0.98 

50,000 

401 

9,892 

0.82 

100,000 

726 

17,921 

0.47 


It should be noted that, the term in £2 containing k w is one or two orders of magnitude smaller than 
that containing hj . Therefore, hj is the dominating parameter. 


COMPARISON AND DISCUSSION OF THE RESULTS 


In order to set forth a precise comparison of the heat transfer predicted by the present model and the 
Nusselt model, one may define two new heat transfer rates with the aid of equations (33) and (35) by 
denoting. 


®Nu = QNu 


f * \ 

u PPg 

-l 

f3nk(T sat -T b )r 0 ) 

“fg 

v ■ ) 


^ ppghfg j 


(41) 


and, 



f \ 

-1/ 

0 = Q 

h f „PPg 

l g 3h J 

V 


3^ik(T sa t Tfr) 

PPghfg 



(42) 


These heat-transfer rates are dimensionless and are closely related to the liquid film thickness at the 
base of the cylinder. As described in the numerical solution section, integration of equation (20), 
provided the values of A for parametric values of £2 . As noted earlier, for £2 = 0 the Nusselt model 
is reached with A^ u = 2.53 113. Using this value, equation (32) can be employed to obtain the 
corresponding values of % , for parametric values of £2 . A Summary of the results is presented in 
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Table 2 for eight parametric values of Q from zero to infinity. The corresponding values of % are 
also listed in the table. 
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Table 2. 


a 

X 

0 Nu 

© 

% difference 

0.0 

1.000 

2.53113 

2.53113 

0.00 

0.5 

0.693 

1.92249 

1.91545 

0.36 

1.0 

0.515 

1.53876 

1.52224 

1.07 

1.1 

0.486 

1.47330 

1.46248 

0.74 

2.0 

0.320 

1.07960 

1.06326 

1.26 

3.0 

0.221 

0.81585 

0.80991 

0.73 

4.0 

0.165 

0.65528 

0.65135 

0.60 

oo 

0.000 

0.00000 

0.00000 

0.00 


It is apparent that the Nusselt model is in notable agreement with the present analysis. As may be 
expected, deviations from the Nusselt model, however small, do appear to exist. The maximum 
deviations are of order one percent and correspond to the moderate values of Q in Table 2. This is 
not a fortuitous outcome. The range of values for various parameters were chosen to span two 
limiting cases. First, as the values of Q. approach zero, the deviations in heat transfer also tend to 
vanish. This may be attributed to progressively larger values of hj which tend to dictate T w to 
approach T b , a constant. On the other hand, at large values Q, h; tends to be small and low heat- 
transfer rates are expected. The wall temperature T w again will assume a constant value, close to 
that of the saturated steam. These limits, of course, can be reached by the influence of all variables 
contained in Q and not just hj . In both limiting cases the cylinder wall temperature assumes a 
uniform distribution and the Nusselt model is approached. The precision by which the Nusselt model 
approximates the present analysis, clearly testifies to the insight often set forth by the old masters. 
Much can be learned by careful analysis of their assumptions and methods. 
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